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Abstract 

Recent experimental results obtained using f3 decay and isomer spectroscopy indi- 
cate an unusual behaviour of the energies of the first excited 2 + states in neutron- 
rich Cd isotopes approaching the N=82 shell closure. To explain the unexpected 
trend, changes of the nuclear structure far-off stability have been suggested, namely 
a quenching of the N=82 shell gap already in 130 Cd, only two proton holes away 
from doubly magic 132 Sn. We study the behaviour of the 2 + energies in the Cd iso- 
topes from N=50 to N=82, i.e. across the entire span of a major neutron shell using 
modern beyond mean field techniques and the Gogny force. We demonstrate that 
the observed low 2 + excitation energy in 128 Cd close to the N=82 shell closure is 
a consequence of the doubly magic character of this nucleus for oblate deformation 
favoring thereby prolate configurations rather than spherical ones. 
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Recent advances in experimental techniques allow to follow the isotopic chains 
of heavy nuclei and therefore the evolution of nuclear structure over wider and 
wider isospin ranges, in some cases even across the entire span of a major neu- 
tron shell. The Cd isotopes (Z=48) for example are the first case, where experi- 
mental information about excited states is now available from the N=50 to the 
N=82 shell closure, that is from 98 Cd to 130 Cd (1; 2). Approaching the proton 
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and neutron drip lines new facets of nuclear structure are expected (3) to gain 
relevance due to the increased importance of the coupling to the unbound con- 
tinuum and changes in the density distributions of protons and neutrons. The 
availability of experimental information over a wide isospin range therefore 
allows for the most stringent test of the validity of global nuclear structure 
models. Besides the nuclear structure aspect, the understanding of the nuclei 
around the doubly-magic tin isotopes 100 Sn and 132 Sn is also of astrophysical 
interest. The rapid proton (rp) capture process of hydrogen burning on the 
surface of an accreting neutron star ends in a closed SnSbTe cycle near 100 Sn 
(4) and the dynamics of the r process and with it the solar abundance dis- 
tribution in the mass 130 region is determined by the properties of the N=82 
waiting point nuclei below 132 Sn (5). Already more than ten years ago it was 
shown that the assumption of a quenching of the N=82 neutron shell closure 
leads to a considerable improvement in the global abundance fit in r process 
calculations (6; 7). Different theoretical calculations indeed predict such a re- 
duction of the shell gap near the neutron drip line (8; 9). Unfortunately, this 
region of very neutron-rich waiting point nuclei is still out of reach experi- 
mentally. However, in recent years a number of experimental observations in 
the 132 Sn region have been interpreted as first experimental evidence of the 
quenching of the N=82 shell closure in 130 Cd (10; 11), much closer to 132 Sn 
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Fig. 1. Excitation energies of the first excited 2 + states as a function of the neutron 
number N for even-even Cd, Te, Xe, Ba, Ce and Nd isotopes in the range N=70-82 
shell. 
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than predicted by any calculation. One of them concerns the anomalous be- 
haviour of the energies of the first excited 2 + states in the heavy Cd isotopes 
as illustrated in Fig. 1. The 2 + state in the N=80 isotope 128 Cd at an exci- 
tation energy of 645 keV lies 7 keV lower than the 2 + level in 126 Cd (10). It 
is evident from Fig. 1 that the neutron-rich Cd isotopes display a completely 
different behaviour when approaching N=82 than the Te, Xe, Ba, Ce and Nd 
isotopes above the Z=50 shell closure. However, the recently observed excita- 
tion scheme of 130 Cd and its comparison to shell model calculations does not 
give any evidence for a reduced N=82 gap in this nucleus (2). Furthermore, 
a new direct mass measurement for 134 Sn restored the 132 Sn shell closure to 
that expected for a doubly-magic nucleus (12). Previous Qp measurements 
seemed to indicate a reduced strength of that shell closure. In light of these 
new experimental results the interpretation of the anomalous behaviour of the 
2 + energies as hint to a possible N=82 shell quenching has to be questioned. 

In this Letter we study the evolution of 2 + excitation energies in the Cd 
isotopes between N=50 and N=82 from a theoretical view point in an attempt 
to unveil the origin of the anomalous behaviour. This endeavor of a realistic 
calculation across an entire major shell can only be realized in the framework of 
modern beyond mean field theories with effective interactions (13; 14; 15). An 
exploratory study of the neutron-rich Cd isotopes has already been presented 
in (16). In that reference, a simplified approach, neglecting some relevant 
exchange terms of the Gogny force, has been used and only a few isotopes 
have been investigated. Furthermore, no particle number projection has been 
performed, neither at the mean field level nor in the configuration mixing 
approach. 

In this work we use the recently proposed symmetry conserving configuration 
mixing approach (17) with particle number and angular momentum projection 
and the finite range density dependent Gogny interaction (D1S parametriza- 
tion (18)) with all exchange terms included. This interaction has been adjusted 
to reproduce global properties of nuclear matter more than twenty five years 
ago and does not invoke any parameter adjustment to certain regions of the 
nuclear chart. Its strength is the ability to describe in a consistent way a large 
variety of phenomena such as the erosion of the N=20 shell closure around 
32 Mg (14), the emergence of the new shell closure at N=32 in neutron-rich 
Ca, Ti, and Cr isotopes (17), the shape coexistence in the Pb isotopes (19) as 
well as the fission barriers in 254 No (20). The intention of these calculations 
is to describe general trends and evolutions of nuclear properties and to help 
to unravel the basic underlying origins rather than to obtain an as detailed as 
possible description of a single nucleus of interest. 

In our approach (17), we first generate the collective subspace of Hartree-Fock- 
Bogoliubov (HFB) wave functions, {\(fi(q))}, by the quadrupole constrained 
particle number projection (PNP) before the variation (21), i.e., we minimize 
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Fig. 2. B{E2) transition probabilities (top) and 2 + excitation energies (bottom) for 
the Cd isotopes between N=50 and N=82. 



the energy 



F N,z (n) = (^ z (g)\H-X q Q\^(q)) 

W (<Z> N ' Z (q)\<Z> N ' Z (q)) ' 1 ' 

with \$ N ' z (q)) = P N P z \(j)(q)), where P N W is the projector onto neutron 
(proton) number and \ q the Lagrange multiplier determined by the constraint 
{tp\Q\tp) / (tp\ip) = q, in an obvious notation. The constraint on the opera- 
tor Q allows the definition of potential energy surfaces along the most rele- 
vant degrees of freedom. In the present calculations we restrict this set to the 
quadrupole deformation (3, triaxial shapes, i.e. the 7 degree of freedom, are 
not considered. 
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The configuration mixing calculation is performed within the generator coor- 
dinate method (GCM) framework taking linear combinations of the particle 
number projected wave functions obtained in the first step and performing in 
addition angular momentum projection: 

1*2?) = / fjfiQ) p j l** z (*)>dg. (2) 

with P J the angular momentum projector. 

Then, the variational principle applied to the weights fjf(q) gives the gen- 
eralized eigenvalue problem (Hill- Wheeler equation): 

/ g') - E»/Af?> z ( q , g')) tf/Mdj = (3) 

with Ti,^' 2 and M^ ,Z the Hamiltonian and norm overlaps, respectively (see 
(14) for further details). In the calculations we use 11 harmonic oscillator 
major shells and the D1S parametrization (18) of the Gogny force. 

It is important to emphasize here that nuclei close to shell closures belong to 
the weak pairing regime where the HFB approach breaks down. To properly 
deal with the pairing correlations in nuclei such as 132 Te and 128 Cd with only 
a few particles/holes outside a doubly-magic core it is mandatory to apply the 
so-called Projection Before Variation approach as described above. As a matter 
of fact we find quantitative differences between the present calculations and 
those of Ref. (16) performed also with the Gogny interaction but without PNP 
and without exchange terms. Furthermore, in order to check the validity of the 
axial approach used in this work, we have also performed triaxial calculations 
at the PNP level. The results show that none of the analyzed isotopes exhibits 
triaxial minima or large 7 softness (22). 

In Fig. 2 we compare the calculated excitation energies of the 2f states (the 
eigenvalues of Eq. 3 for J = 2 + and o = 1) and the transition strengths 
B(E2, 0+ -> 2+) for all Cadmium isotopes from the N=50 to the N=82 shell 
closure to the available experimental information (23). We find an overall good 
qualitative agreement between theory and experiment. The general trend as 
well as the absolute strength of the experimental B(E2) values are well repro- 
duced by our calculations. For the 2 + excitation energies one can distinguish 
three different regions. In the first one, N=51-58, the parabolic behaviour 
of the energy roughly corresponds to the successive filling of the gr/2 neu- 
tron orbital. In the second region, N=59-70, with a rather flat energy, several 
smaller sub-shells (d 5 / 2 , d 5 / 2 , Si/ 2 ) are being filled. Finally, from N=71 to 82, 
the hn/2 orbital is being filled, but at variance from the first region here we 
observe the unexpected non-parabolic behaviour of the 2 + energies. These 
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Fig. 3. a),b) Calculated and experimental 2± energies in the Cd and Te isotopes 
as a function of the neutron number; c),d) potential energy curves as a function 
of the intrinsic quadrupole deformation f3, with particle number projection (dotted 
lines) and after additional angular momentum projection, J = Oh (continuous lines) 
and J = 2h (dashed lines); e),f) squared amplitudes of the wave functions of the 
collective states of Eq. 2 with J = Oh (continuous lines) and J = 2h (dashed lines). 

three regions have an approximate counterpart in the transition probabilities. 
Though we nicely reproduce these features the calculated energies are too high 
on an absolute scale. The origin of this scaling is supposedly due to the ab- 
sence of K- mixing in the calculations. This mixing can be induced in different, 
non-independent, ways: at the HFB level, by angular momentum projection 
before variation (24), by the Coriolis force (cranking) or by admixing of two- 
quasiparticle states (Tamm-Dancoff) (24; 25), and at the configuration mixing 
level by the consideration of triaxial shapes (26). Within our approach we have 
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investigated the (ft, 7) plane and, as mentioned above, none of the nuclei stud- 
ied show triaxial minima or large 7 softness and it seems that in this case the 
consideration of K-mixing effects amounts, at least for the lowest lying collec- 
tive states, to a global, rather constant, change in the mass parameter (either 
for rotations or vibrations). That means, the ratio between the calculated and 
experimental 2 + energies is nearly constant (about 1.6) over a large range of 
isotopes. This effect, which is not yet fully understood, has already been ob- 
served in previous studies, see for example (17). A rather constant ratio has 
also been obtained in limiting situations evaluating the mass parameter with 
and without triaxial effects (27) and in an exact triaxial angular momentum 
projection (28). Taking into account this constant factor the calculations re- 
produce the experimental data even quantitatively. This is shown in Fig. 3(a,b) 
where the renormalized calculated 2 + excitation energies for the Cd and Te 
isotopes are compared to experimental data in the region of interest close to 
N=82. The calculations nicely reproduce the experimental trends, namely the 
parabola-like behaviour in the Te isotopes and the "pathological" flattening 
up to N=80 for the Cd isotopes. 

To trace back the origin for the different behaviour of the N=76, 78, and 80 Cd 
and Te isotones we will compare in more detail the different steps of the calcu- 
lations for the N=80 pair. In Fig. 3(c,d) we present the potential energy curves 
calculated before configuration mixing as a function of the intrinsic quadrupole 
deformation ft. These curves show clear differences between 128 Cd and 132 Te 
and shed light on the origin of the different behaviour of the 2 + energies. The 
energies obtained using the solutions of the PNP equations (see eq. 1) show an 
intrinsic spherical minimum in 132 Te and a slightly prolate one in 128 Cd. The 
results obtained after angular momentum projection, i.e., the energy expec- 
tation values calculated with the wave functions Z (q)) = P J P N P z \(j)(q)), 
are rather similar for both nuclei for the + state: They present coexisting 
prolate and oblate minima. However, the results for the 2 + state are rather 
different: for 128 Cd we obtain a deep prolate deformed minimum about 3.5 
MeV lower than the oblate one, whereas in 132 Te we observe again two sym- 
metric minima as for the + state. Finally, in Fig. 3(e,f) the wave functions of 
the collective states (eq. 2) with angular momentum and 2, respectively, i.e. 
the states obtained after configuration mixing, are shown (the corresponding 
2 + energies have already been shown in Fig. 3(a,b)). Once more, pronounced 
differences are observed for the 2 + state: In 132 Te this state represents on 
average a quadrupole oscillation around a spherical shape. In 128 Cd, on the 
other hand, the 2 + state is essentially an oscillation around a prolate deformed 
shape and therefore strongly coupled to rotation. This is the main reason why 
the 2 + energies in the heavy Cd isotopes behave differently from those in 
the heavy Te isotopes. Trying to understand the underlying origin of the dif- 
ferences between the properties of the 2 + states in 128 Cd and 132 Te we will 
now have a closer look to the potential energy curves obtained after particle 
number but before angular momentum projection (dotted lines in Fig.3(c,d)). 
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Fig. 4. Proton (left) and neutron (right) single particle energies for 128 Cd. The 
energy scale for neutrons has been shifted by 9 MeV. The thick dotted-dashed lines 
represent the Fermi levels for 128 Cd and the thick dashed line the proton Fermi level 
for 132 Te. 

Already here, a difference between the two nuclei is apparent. On the oblate 
side {(3 < 0), the energy increases much faster for 128 Cd compared to 132 Te 
(E(/3 = -0.1) - E((3 = 0)=2.5 MeV and 1.0 MeV, respectively). This dif- 
ference can be understood on the basis of the single particle energies (SPE), 
represented in Fig. 4 for 128 Cd. The SPE for 132 Te are not significantly differ- 
ent from these ones. The proton Fermi levels for both nuclei are included in 
Fig. 4, while for the neutrons we assume the Fermi level to be the same since 
both nuclei have the same number of neutrons, N = 80. At an oblate defor- 
mation of about (3 = —0.1 the neutron lsj/2 level is emerging from the Fermi 
surface, see the corresponding arrow. From here on both N=80 isotones have 
a completely filled z^/iii/2 shell. On the proton side, a similar situation occurs 
around (3 = —0.15 where the proton lpi/2 level is crossing the Fermi surface 
for 128 Cd. That means that for oblate deformations larger than \f3\ m 0.15 in 
addition to the full hn/2 neutron shell the nucleus 128 Cd also has a completely 
filled large intruder proton orbit, namely the 7rg 9 /2 shell. To deform this nu- 
cleus with closed large intruder shells for both protons and neutrons a lot of 
energy is required resulting in the very steep PNP curve for 128 Cd in Fig. 3(c). 
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Furthermore, the generation of two units of angular momentum to form a 2 + 
state is very costly for such a stiff configuration since it requires the admixing 
of orbitals above the Z = 50 or iV = 82 shell closures (compare the curve for 
J = 2 + curve in Fig. 3(c)). The situation in 132 Te is very different. In this 
case, the proton Fermi level lies above the Z = 50 shell closure and the oblate 
branch of the proton system can easily be deformed (see Fig. 4) due to the 
positive slope of the relevant #7/2 orbitals. This deformation will be directly 
transfered to the neutron system due to the strong proton-neutron interaction. 
Even more important, a J = 2 + state can easily be build by coupling two #7/2 
protons (see Fig. 3(d)). To summarize this discussion, we have shown that the 
prolate deformation of the 2 + state in 128 Cd has its origin in a very peculiar 
coincidence of closed shell proton and neutron configurations for {3 < —0.15 
leading to a "blocking" of the oblate branch of the energy curve, i.e., in this 
branch the nucleus behaves as a double magic one. Since the proton oblate 
blocking persist also in the nearby iV = 76 and 78 Cadmium isotopes it is 
obvious that these isotopes will be on average more deformed than the cor- 
responding Tellurium isotopes, as depicted by the 2 + energies in Fig.3(a,b), 
indicating that the different behavior of the Cadmium and Tellurium isotopes 
close to the iV = 82 shell closure can be explained by "standard" nuclear 
structure effects without invoking shell quenching. 

In conclusion we have presented a theoretical study of the origin of the unusual 
behaviour of 2 + excitation energies in neutron-rich Cadmium isotopes towards 
the N=82 shell closure using modern beyond mean field techniques and the 
Gogny force. We have demonstrated that these calculations, which do not 
dispose of any adjustable parameters, are not only able to describe general 
features of atomic nuclei across entire major shells, i.e. over wide ranges of 
isospin, but also to reproduce local features of nuclear structure such as the low 
2 + energies in 126 > 128 Cd and unravel their origin. Furthermore, our calculations 
show that the anomalous behaviour of the 2 + energies in the heavy Cd isotopes 
is caused by the very special characteristics of the Cd isotopes which favor 
prolate configurations close to the N = 82 shell closure. There is no need to 
assume any changes of nuclear structure far-off stability, i.e. a quenching of 
the N = 82 shell closure, to explain the experimental findings. 
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